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Abstract 

' A statistical theory of the mean field is developed. It is based on the propo- 

' sition that the mean field can be obtained as an energy average. Moreover, 

it is assumed that the matrix elements of the residual interaction, obtained 
^ , after the average interaction is removed, are random with the average value 

\ of zero. With these two assumptions one obtains explicit expressions for the 

<0 ' mean field and the fluctuation away from the average. The fluctuation is 

expanded in terms of more and more complex excitations. Using the random- 



o 



■ ness of the matrix elements one can then obtain formulas for the contribution 

to the error from each class of complex excitations and a general condition 
for the convergence of the expansion is derived. It is to be emphasized that 
no conditions on the nature of the system being studied are made. Making 
some simplifying assumptions a schematic model is developed. This model 
^ ■ is applied to the problem of nuclear matter. The model yields a measure of 

the strength of the effective interaction. It turns out to be three orders of 
^ \ magnitude less than that calculated using a potential which gives a binding 

5h \ energy of about — 7MeV/nucleon demonstrating the strong damping of the 

interaction strength induced by the averaging process. 
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I. INTRODUCTION 



Remarkably, the mean field provides a good description of many systems. The underlying 
dynamics may be very different, the interparticle interaction may be long range or short 
range; it may be weak or strong. Nevertheless in each case a mean field description agrees 
closely with experiment. The question naturally arises: what are the common features which 
these various systems have which ensure the effectiveness of the mean field approximation? 
This is the issue which mostly concerns us in this paper. 

We proceed by presenting a formalism in which the mean field and the error are derived 
without any reference to a particular dynamical system. The mean field is obtained as 
an energy average which smooths over short time events. This procedure also develops an 
expression for the error which must be small if the mean field is to be useful. To calculate 
the error a related assumption is used. Not only is an energy average taken, but it is also 
assumed that the error matrix elements are random with the consequence that the average 
error is zero, but of course the average of the square of the error is not. The error can be 
calculated as a finite series consisting of contributions from various degrees of excitation. 
A general expression for each term is developed and a parameter determining the relative 
importance of each term is identified. 

It will be noted that the assumptions listed above, namely energy averaging and random 
error matrix elements, are identical with those used in the statistical theory of nuclear 
reactions [Q. 

This paper concludes with an exploratory calculation appropriate for the ground state of 
nuclear matter. Only the first order error is considered. One obtains reasonable estimates 
of the error. But importantly there is an enormous reduction in the strength of the residual 
interaction from that obtained from a Fermi gas model which must be a consequence of the 
statistical approach. 

The statistical approach provides a method for calculating the slowly varying and rapidly 
varying contributions to physical observables. It differs from the customary approach based 
on the Pauli principle and dynamical correlations. The latter are a consequence of the 
rapidly varying contributions. The theory presented in this paper is essentially a different 
way of taking these contributions into account. 

Our work is organized as follows. In sect. |T| the formalism is revisited, — in particular 
the partition of the Hilbert space into the P and Q sectors, — and the related equations 
are discussed, while in sect. Ill a precise definition of the energy average is provided. The 
latter is performed with a distribution, characterized by a single parameter e, amounting to 
shift the position of the poles corresponding to the excited states of the system away from 
the ground state energy. In the same section, we pave the way for setting up a suitable 
formalism to account for the corrections (the error) to the mean field energy. 

In sect. fV| the expression for the corrections to the mean field energy Eq is given in the 
framework of the expansion, above referred to, built out of successive 2p-2h excitations. 

In sect. |V] we deal with the issue of introducing an operator, projecting into a smooth 
P-space, simple enough to allow the present approach to be worked out for nuclear matter. 
We thus restrict the P-space, — also with the aim of rendering transparent the comparison 
between the mean field energy of the Hartree-Fock (HF) and of the present theory, — just 
to a simple state: The HF determinant. With this choice for P, we are able to establish 
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an expression connecting the true ground state E, the mean field Eq and the Hartree-Fock 
i?HF energies. Importantly, Eq is shown to be always lower than E^p and it turns out to 
be proportional to the variance, to be later defined, of the residual effective interaction V 
among nucleons, which, in our context, is precisely defined. 

In sect. |V| we address the issue of the occupation probability S"^ (spectroscopic factor) 
of the ground state of nuclear matter as defined by the structure of our P-space. We are 
able to deduce, within our scheme, an expression for S*^ in terms of E, Eq and -Ehf, and 
moreover embodying the energy derivative of V . 

Finally, in sect. |Vli| a scheme which permits numerical predictions of the present approach 
is developed. 

For sake both of simplicity and of illustration, we account for the first corrections only 
and, even so, we are forced to introduce a free parameter a in order to give an estimate of 
their size. We are then able to set up two systems of two equations each by coupling the 
equation for the mean field energy to the one for the fiuctuations. By solving the systems 
and with a convenient choice of the two parameters entering into our approach, namely e 
and a, we obtain results consistent with the empirical features of nuclear matter as obtained 
through an extrapolation from finite nuclei. We provide at the same time an estimate for 
the fiuctuations of the mean field energy. These goals are achieved with a residual effective 
interaction which, as previously emphasized, is lowered by about three orders of magnitude 
with respect to the bare one and yields a spectroscopic factor not in conflict with other 
independent estimates 0. 

In the concluding section our results are summarized and possible improvements dis- 
cussed. 



II. FORMALISM 

In this section we shortly revisit the formalism introduced in Ref. 0. Let H be the 
nuclear hamiltonian entering into the Schrodinger equation 

Hm = Em. (2.1) 

Let P be the hermitian operator projecting the nuclear ground state into the Hilbert sub- 
space of functions associated with the low momentum transfer physics and Q the operator 
complementing P. Clearly 

p2 = p, g2 = g, pg = o, p + g = i. (2.2) 

It is then proved that the pair of equations 

{E - Hpp) (P^) = HpQ (g^) , (2.3) 

and 

{E-HQQ){Qm) = HQp{Pm), (2.4) 

are equivalent to (|2.1| ). In the above formulas the following shorthand notations have been 
introduced, viz. 
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Hpp = PHP, Hqq = QHQ, HpQ = PHQ, Hqp = QHP. 



(2.5) 



Let us next define an auxiliary function \l/o, wliicli in tlie end will disappear from the 
formalism, such that 



{E - Hpp) ^0 = 0. 



Then, from (|2.3|) it follows 

(P^) = ^0 
which, when inserted into 

(W = 

being 



E-H. 



pp 



Hpq{Q^) 



yields 
1 



E-Hoo- W, 



Hnp^O = Hnp \&o 



QQ - y^QQ 



Wqq = Hqp 



E — Hpp 



H 



PQ 



and, obviously. 



E - Hqq - Wqq . 



By using ( p.7| ) and ( |2.8| ) one can, on the one hand, recast ( p.3|) as follows 



{E - Hpp) (P^) = HpQ — Hqp *c 



and, on the other, express \&o according to 

1 



1 



(Pvl/). 



E — Hpp 



HpQ — Hqp 



The insertion of ( |2.12| ) into leads in turn to the equation 



{E - Hpp) (P^) = HpQ — Hqp 

eg 



which can be reorganized into the form 

1 



1 



1 + FT — HpQ — Hqp 

h, — Jrlpp Cq 



{E - Hpp) (P^) = HpQ 



1 



1 + Hqp ^ _ ^ — HpQ 



Hqp{P^) 



HpQ — 



pp 

Hqp (Pvl>) . 



^ + Wqq — 



(2.6) 



(2.7) 



(2.8) 



(2.9) 



(2.10) 



(2.11) 



(2.12) 



(2.13) 



(2.14) 
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In deriving ( |2.14| ) the operator identity 

1 



B 



1 



1 + CAB 1 + BCA 



B 



has been used, setting 



A 



B = H, 



QP ; 



C 



E — Hpp 



H 



PQ 



Then, since 



one finally obtains the equation 
/ 

E — Hpp — H 



eg. 



QQ 



PQ 



V 



— +Wc 



QP 



(P^) = 0. 



We thus see that the intricate many-body operator 

1 



Tt = Hpp + HpQ 



H 



eg. 



QP J 



QQ 



(2.15) 



(2.16) 



(2.17) 



(2.18) 



(2.19) 



dependent upon the exact ground state energy of the system E, defines an equation which is 
not an eigenvalue equation in the usual sense. Indeed, the dependence upon E is non-linear, 
since the latter appears in the propagator in the Q-spa.ce in (|2.18| ). 



III. ENERGY AVERAGE 

As emphasized in the introduction, {P"^) contains both slow and rapidly varying com- 
ponents. To obtain the slowly varying components one must perform an energy average of 
the wave function. The components varying most rapidly with energy are taken to be in the 
Q space so that {Q"^) is to be replaced by {Q"^) . 

As a consequence of this procedure also the {P"^) component of the wave function and 
the energy E as well, will no longer be the same (see eq. ( p.3| )). We shall call the modified 
quantities {P"^) and Eo respectively and obtain, instead of (|2.7|) , the following equation 

(Pvl/) = ^0 + ^ HpQ {Q^) , (3.1) 

tLo — tipp 

where 
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[Eo - Hpp) ^0 = 0. 



(3.2) 



Analogously, eq. ( p.8|) will now read 



(Q^) = (— ) Hqp^o 

eg 



which, upon insertion into ( p.l|) , leads to 



Eo-Hpp]{P^) =Hpq{ — ) Hqp^o. 



(3.3) 



(3.4) 



Exactly as before, '^o can now be obtained by combining (|3.1| ) with ( p.3|) and, when inserted 
into (|i|), it yields 



/ 



pp 



H 



PQ- 



-1 



QP 



\ 



(— ) + 



QQ 



(P^) = (Eo - h) (P^) = 



(3.5) 



which is the correspondent of equation (|2.18| ) when a suitable energy average has been 
performed. 

At this point we shall specify how the energy average is to be carried out. One introduces 
a smoothing function p{E, Eq) with the property 

(3.6) 



(3.7) 



p{E, Eq) dE = l. 
The average of a function f{E) is then 

(/) = I p{E,Eo)f{E)dE. 



Since we are considering bound states we require the average (/) to be real. A smoothing 
function which satisfies these conditions is 



p{E,Eo 



1 



1 



27ri E-{Eo- e) 



(3.8) 



with a path of integration in eq. (|3.7| ) going along the real axis ReE with a small semi-circle 
described positively about the singularity (^o ~ Assuming boundedness conditions at 
infinity sufficiently strong, the Cauchy's integral formula can be applied so that the condition 
of eq. p.6| ) is satisfied and eq. p.7| ) becomes 

(/) = f{Eo - e) . (3.9) 

Following the procedure developed by Kawai, Kerman and McVoy 0], we find 



(p - n) (p^) = VpQ (Q^) 

{E-HQQ)iQ^) = VQP (Pvj/) 



(3.10) 
(3.11) 
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where 



Ti = Hpp + VpqVqp 



Eq — e — E 



(3.12) 



and 



V, 



PQ 



H 



PQ 



Ep-e-E 



QP 



E 



He 



QP ■ 



(3.13) 
(3.14) 



Note that eqs. ( |3.1CI| ) and ( |3.11| ) are of the same form as the original eqs. ( |2.3| ) and 
However the potential which couples {P'^) and {Q'^) is Vpq rather than Hpq. The strength 
of the coupling is thus considerably reduced by roughly (e/HQq)^ which is much less than 
one. This assumes that Hpp and Ti, are of the same order of magnitude. 

Concerning the parameter e, in principle its value should be inferred from experiment; 
in practice we shall consider it as a free parameter to be utilized in the search for agreement 
with experiment. 

Now the spectral decomposition of the operator 1/{E — Ti) can be performed: 



1 



E-n 



( $0 I VpQ I ) 
E — Eq 



E — Er, 



+ 



E-n 



Vpq{Q^) 
Vpq{Q^) 



(3.15) 



where 



n — En^n 



(3.16) 



The prime on {1/{E — Ti.)) signify that the lowest eigenfunction $o is to be excluded. 
From eq. ( |3.15| ) we have 



( $0 I ) 



$0 I VpQ I ) 



E — Eq 

To obtain {Q"^) we return to eqs. ( |3.10| ) and ( p.ll[ ) and get 

1 



E-h 



VQp\%){%\P<i/), 



where 



and 



QQ 



Hno + W, 



QQ 



QQ 



Vc 



QP 



1 



E-n 



PQ 



(3.17) 



(3.18) 



(3.19) 



(3.20) 
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g. s. 2p-2h 4p-4h 

FIG. 1. The partition of the Hilbert space of nuclear matter in sets of increasing complexity. 
The first box on the left defines the P-space, the second one embodies the simplest states in the 
Q-space and so on. 



Inserting eq. (|3.18|) into eq. (|3.17| ) one obtains 



E — Eq = { \ VpQ 



1 



E-hc 



QP 



$0 ). 



(3.21) 



In this equation Eq is the mean field energy while the right hand side represents the correc- 
tion. 

Note that the quantity (•I'ol-P^) does not enter into eq. ( |3.21| ). This is to be expected 
since, so far, the amplitude of $o is arbitrary. However since $0 is the wave function for a 
bound state it can be normalized. In the present context we choose ($o|-P^) = 1- 



IV. CORRECTIONS TO THE MEAN FIELD ENERGY 

The corrections to the shell model energy Eq are obtained in the statistical formalism 
as in the case of the nuclear reactions by expanding ( |3.21D in a series of contributions 



which, rather than by powers of a coupling constant, are identified by the complexity of the 
states of the Q-space. 
Thus one writes 

AE=j2^Em, (4.1) 



m=l 



where 

^lol^o), (4.2) 
= Q,VQ, , (4.3) 

and Qj is the operator projecting into the j particles-j holes sector of the Q-space. 

One thus sees that the states belonging to the Q-space are classified into r sets, each set 
incorporating a well-defined type of excitation; specifically the k-th set embodies all the k 
particles-/^ holes states. Actually, because of the essentially two-body nature of the residual 
interaction V, the one effectively responsible for the passage from one set of states to the 
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other, only excitations corresponding to an even number of particles (and holes) are actually 
filling the sets, as indicated in Fig. |l|. The propagation of the system into the k-th set is 
then described by the operator Gk, which when k = r reads 

Gr = y^ , (4.4) 
whereas for k < r obeys the recurrence relation 

E — hkk — Vk^k+iGk+iVk+i^k ^ ^ 

In the above 



hkk = Hq^q^ + Vq^p i^^ _ VpQ^ . (4.6) 

For the sake of illustration each set is sometimes displayed as a box: The order among the 
sets, which is a condition for the validity of (|4.2|) , can then be expressed by stating that no 
box can be occupied unless the previous one has already been too (see Fig. p. 

Moreover, it should be noticed that in a nucleus the expansion ( [4.1|) is finite since the 
excited states with the greatest complexity correspond to the case r = A, A being the mass 
number. 

Since the mean field contains the average, the average of the corrections to the mean field 
vanishes. In addition because of the random nature of the processes occurring beyond the 
mean field the average of the corrections to the shell model stemming from any individual 
set vanishes as well. Thus 

{AE) = {AEi) + {AE2) + ■■■ + (AEa) = (4.7) 

and 

(AEn) = 0. (4.8) 

However the average of the square of the corrections does not vanish. For this, namely for 

{{AEf) = {Y,AEnAEm), (4.9) 

n,m 

the randomness of the Q-space physics is exhibited in the randomness of the nuclear matrix 
elements entering into the definition of AEn- As a consequence, the off-diagonal elements 
in (|4.9| ) disappear, leaving 

{{AEf) = {Y.{^E^?). (4.10) 

n 

The implementation of the requirement 

(A^i) = {AE2) = ■■■ = (AEa) = (4.11) 
is achieved, in the case r = 1, by redefining AEi itself according to 
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AE, = ( 00 I K)iGiVio I 0^ ) - [( 0^ I VoiG^V.o \ 0^ ) 



AV 



(4.12) 



and likewise for the terms with r 7^ 1. In ( |4.12| ) 0^4 represents the nuclear mean field ground 
state wave function and the square brackets mean energy averaging. 

Now, the corrections to the mean field can be calculated, according to the formula ( p. 211 ), 
provided the eigenfunctions tp^^a and the eigenvalues ek^a of the operator are known, 
which in turns implies that the equation 

Gk^i'k,a = {E - ek,a)^k,a 

has to be solved for any value of the index k. 
To begin with we consider AEi: 



(4.13) 



f3 E-e2i3 



E 

/3 



E — €2(3 



(4.14) 



AV 



where €2/3 and |'02/3) are the eigenvalues and eigenvectors, respectively, corresponding to the 
2p-2h excitations. 

However, as previously remarked, in the present scheme only the square of the corrections 
to the mean field, namely 



(5^1 



AE,f - [(AEi)] 



AV 1 



(4.15) 



are meaningful. We are left with the task of evaluating 



{5E, 



E 

/37 



AV 



{E - t2p){E - €2.) 

For the first term of the above equation we have 

(<^^|V^|^2a)(^2a|l^|0V(0Al^l^2/3)(^2;3|^|0V 



fi E - €2/3 



(4.16) 



AV 



E 

Q!,/3 



{E - e2o){E - e2fi) 



AV 



E-t2 



E(</'a|V^I^2«)(^2.|^|0^)(0^|V^|^2/3)(^2/3|1^|0^) 

a,/3 



(4.17) 



AV 



where €2 is the average excitation for the 2p-2h states. To obtain the random phase average 
of the quantity in the square brackets we make use of the result: 



(v4v4*S5*) = {AA*){BB*) + {A*B){AB*) + quartic terms. 



(4.18) 



We shall neglect quartic terms. 

The first term {^AA*){BB*) cancels exactly the second term in eq. ( [4.161) . Evaluating 
the second term in eq. ( [4.18|) yields 
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2a/ Uy 



2 

AV 



(4.19) 



In this last equation AEi is the interval over which the energy average is carried out. It 
is essentially equal to e. The quantity D2 is the energy distance between two neighboring 
2p-2h states. Finally ip2 is meant as a representative of the 2p-2h excitations. 
We now go on to consider AE2: 



AE2 = (0^1 VoiG'lVl2G2y2lGi\/io 10^) 



E 



E — e2a E — 64/3 
1 1 



E — e2a' 



{E-e2)HE-e^) 

2a I 



E(</'A|V^IV^2a)(V^2.|V^|0^) 



a,l3 



AV. 



(4.20) 



where we have set a = a' and in the last line have subtracted the average of AE2 as written 
in the first line so that {AE2) = 0. We now consider the square of the modified AE2 = SE2. 
We obtain 



{6E2 



1 



E {<pA\V\^p2a){^2a\V\<P'^) 



(E - 62)4 (E - 64)2 

(V^2a|V^|^4/3)(V^4/3|V^|V^2a)^0^|V^|^2a')(^2a'|^|0^) 



{^2a' I V I V'4/3' ) i^Pifj' I V I ?/'2a' ) 
|(0°|^|^2a')ri(^2a|^|^4/3)ri(^2a'|V^|^4;3')l 



E [\{<pA\v\t 



'2a 



AV 



AV 



(4.21) 



To calculate the average we now invoke eq. ( |4.18| ). The first term again cancels the subtracted 
term and one is left with 



(5^2 



(E - -e2Y (E - 



EK</'^l^l^2.)n(^2.|V^|^4,)r 

Q,/3 



(4.22) 



AV 



Comparing this result with that for {6Ei)'^, eq. ( |4.19| ), one finds 

1 1 



{SE2)' 



{E - £2)2 {E - €4^2 



{6E,y\{Hv\^,K- 



(4.23) 



which gives us an expression for the expansion parameter, though, as it has been emphasized 
earlier, the expansion has a finite number of terms. Again in the above 1^4^ is portraying 
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a representative of the 4p-4h states. Eq. ( |4.23| ) can be generalized into a relation between 
{SE^f and 



{5E„ 



(E - e^n? (E - -e2n-2? 



{6E^^,f 



|(^2n-2|^|^2n)r 



(4.24) 



This concludes the general development of the statistical theory. To make further 
progress a detailed description of the system under consideration is required. Our appli- 
cation will be to nuclear matter. But one should bear in mind that up to this point our 
results can be applied to any system to validate the mean field concept. 



V. THE PROJECTION OPERATOR AND THE HF THEORY 

To obtain the true energy {E) and the mean field one (^o) (the shell model energy for 
a nucleus) we now need an explicit form for both the projection operator P and the NN 
interaction. 

Let us start first by searching for a suitable projection operator. Suppose, for this 
purpose, that for a given interaction one has solved the associated HF problem. Let us call 
{0^p}, with the index i = 1,2, . . ., the single particle wave functions of the HF orbitals, 
which form an orthonormal complete set for the Hilbert space of a one particle system. Out 
of the {0hf} ^6 can build an infinite set of Slater determinants {xhf} which also form an 
orthonormal complete set for the Hilbert space of the nucleus. 

We can then define 

M 

^ = ElxW)(xkFl (5.1) 

i=l 

which, in the special case of M = 1, reduces to 

P=\xIif){xIif\, (5.2) 

IXhf) being the ground state HF determinant of the nuclear system. In the following we 
shall stick to the option M = 1, although generalizations to larger values of M, while 
cumbersome, should be useful to consider. 

Notice that, with P given by (|5.2| ), the wave functions [P"^) and (P"^) turn out to be 
proportional to each other and to the ground state HF determinant |xhf) they must since 
the P space has only one member. 

With this projection operator and the mean field given by H (eq. (|3.12| )) 



Eo = EnF + { Xhf \VQQV^ Xhf )^ ^ (5.3) 

where 

EnF = ( Xhf I Hpp I Xhf ) . (5.4) 

and 
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^ - e - Hqq 

(see eq. ( p.l3| )). Eq. ( ^.3|) , writing Q = 1 — P, becomes 

^0 = + ^^_\_^ {{ XuF I VV^ I xJiF ) - l( Xhf I V I Xhf )r} • (5.6) 

The expression ( p. 6] ) is general. From now on, however, we shall confine ourselves to 
consider a Fermi liquid (infinite nuclear matter). For such a system |xhf) = 1-^)) 1-^) being 
the wave function of a Fermi sphere of radius kp (the Fermi wavenumber). 

In infinite nuclear matter the quantity in curly brackets in eq. ( |5.6| ), which we refer to 
as the variance of the residual interaction V, can be simplified. One has (the states |n) 
labelling the spectrum of the Fermi sphere) 

{F\ VV^ \F) = V \n){n\ \F) 

n 

^ {F\V\F){F\V^F) 

Ip-lh 

+ {F\V\2p-2h){2p-2h\V^ \F) 

2p-2h 

= \{F\V\F)\' + P\ 

the piece related to the Ip-lh states essentially vanishing as required by the Brillouin theorem 
and further terms in the right hand side being neglected because of the prevalent two- 
body character of the residual interaction V. The quantity (3"^ is just the value of the variance 
given in eq. ( |5.(j| ). The latter is of course hard to deal with because it is not only energy 
dependent, but dependent upon the energy averaging procedure as well. Furthermore, its 
connection to the Q-space Hamiltonian is highly non-linear. However, in nuclear matter 
the variance of the residual interaction is positive since (defined by the eq. ( p.7| ) itself) is 
positive. Moreover, by inserting (|5.71 ) into (|5.6|) , one obtains the important expression 



Fo = FuF + -^r^ , (5.8) 

which embodies all the energies characterizing our problem, namely the HF, the mean field 
(which would correspond to the shell model in a finite nucleus) and the true one. 

Thus, even if a residual interaction V is given and the associated HF energy F^p is 
calculated, still the energy F of the system cannot be obtained since the mean field energy 
Eq remains to be fixed. Yet, if as a first orientation one sets Eq ~ E, then the remarkable 
result 

Eq = Ehf - y (5-9) 

follows, showing that the mean field energy is lower than the HF one, the parameter e being 
of course positive. This result holds only if |^o — <^ e. 
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VI. SPECTROSCOPIC FACTOR 



The projected wave function P\I/ is just a component of the full wave function \E'. Just 
how big a fraction of is contained in P\l/ is obviously of interest. Toward this end let 
(P^|P^) = S"^. We will now obtain an equation for S"^ in terms of the other quantities 
characterizing our problem. We begin with the identity 



= 1 - {Q^m . 



Using eq. ( |3.18| ) for (Q^P), 5^ becomes 

S' = l-{%\VpQ 



Vqp\%) 



We now take the Fermi gas model to be written |P). Then $o = S\F) and 

1 ' 



S' = l-S'{F\VpQ 



Vqp\F) . 



Now, a little algebra allows us to recast the previous equation into the form 



d 



S'^l + S\F\VpQ— 



1 



dE \ E-h 



QQ, 



d 



-{F\ 



dE \ 
dVpQ 



V, 



PQ 



E-hc 



Vqp\F) 



-Vqp \F) 



dE E-h 



^ Vqp\F) - {F\VpQ- ^ 



dVqp 



QQ 



E-hr 



dE 



where the weak energy dependence of hqq has been neglected. Moreover, since 

dVpQ 1 



dE 



-Vj 



PQ 



2{Eo 



E) 



and 



dVr 



QP 



dE 



-Vc 



QP 



2{Eo 



E) 



d 1 



it follows that 
or, by virtue of ( |3.21D , that 

^2 = 1 + ^2 
= 1-S' 



Vqp\F) + 



En 



E 



{F\Vi 



PQ 



E-h 



QQ 



-Vqp\F) 



Je^^ ~ + E^,- e - P 
'dEp 6 \ 

dE Eo-e-Ej ' 



(6.1) 



(6.2) 



(6.3) 



(6.4) 



(6.5a) 



(6.5b) 



(6.6) 



(6.7) 
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The explicit dependence of Eq upon E is provided by ( |5.8| ) which, when inverted, yields 

(6.8) 



En 



E + e + EnF± J{E + e - E_ 



)2 + 4/?2 



where the minus sign in front of the square root should be taken, since for vanishing residual 
interaction (/J^ = 0) the mean field should reduce to the HF field (^o = -^hf)- Hence 



dEo 
dE 



1 lEnF-E -e 



2 dpydE 



2 2 hE^^_E-e)^ + 4p^ 



(6.9) 



which, inserted into (|6.7| ), finally leads to 



3 lEuF- E-e-2 dp^/dE 
2^2 ' ^ 



(Ehf -E-ey + A(3' 



3 1 -Ehf 
2 2.I^E, 



E 



2 _L z1/^2 

2 d[5'^/dE 



Eo-E-e 

Eq — EuF 



HF 



E- 



/32 



(6.10) 



namely to the equation we were looking for. Eq. ( |6.10| ) expresses the spectroscopic factor 
S"^ in terms of the HF, mean field and true energies, with an additional dependence upon 
the energy averaging parameter e, the residual interaction /?2 the derivative of the latter 
with respect to the energy. Since all these quantities are either explicitly evaluated (-Erf, Eq 
and E) or fixed, in principle, both by the experiment (e) and by our theoretical framework 
(/3^), eq. ( |6.10| ) provides an important check for the consistency of our approach, also because 
independent estimates on the depletion of the ground state wave function are available for 
lead, a nucleus to which the nuclear matter concepts should apply. 



VII. A SCHEMATIC MODEL 



In this section we shall employ a schematic model to gain an insight into the effectiveness 
and self-consistency of the formalism described above. The model will rely on two equations 
which, for sake of convenience, are repeated here. The first of these is eq. 



Eq = EuF + 

The second is obtained from Eq. (|4.19|) , i. e. 
which is approximated by 



/3' 



En — E — e 



E-e2 



E[I('^aI^I^2,)I 



2 

AV 



{E - ae)2 



(7.1) 
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where a is a parameter which estimates the average energy of the 2p-2h states involved (we 
set €2 = ae for convenience) and importantly helps to correct for the poor approximation of 

for the sum over the 2p-2h states 7 in eq. ( ^.191 ). Taking the square root of eq. 
yields 

E-Eo = . (7.2) 

E — ae 

If now one combines eq. (|5.8| ) and (|7.2|) , using the upper sign (plus) one obtains a lower 
bound El for the energy of nuclear matter, while using the lower sign (minus) one obtains 
an upper bound Eu- Moreover, one obtains differing values of the mean field energy to be 
denoted by Eq and E^, respectively. We ask whether there is a choice of e, a and j3 such 
that 

i) the two mean field energies Eq and Eq turn out to be the same over a range of densities 
(or of Fermi momenta kp) of significance for nuclear matter; 

ii) the "experimental" values of the binding energy, saturation density and compression 
modulus of nuclear matter, namely B.E./A = —16 MeV, kp = 1.36 fm^^ and K = 23 
or 16 MeV (the two values refer to a hard and a soft equation of state, respectively) 
[0] are accounted for in a sense to be later specified; 

iii) the spectroscopic factor given by eq. ( |6.10|) turns out to be meaningful, in the sense of 
being of course less than one both on the lower bound (where its value Si is associated 
with El) and on the upper one (where its value is associated with Ey_). Furthermore 
Si and Su should not be too much at variance with existing estimates. 

To satisfy condition i) we proceed by first eliminating Eq and Eq from eq. ( p.8| ) and 
This yields an equation for Pf and 



A' = ^^[2El - e{a + 1) - Euf] - ^[2Ei - e{a + 1) - Euf? + MEi - Euf)] 



(7.3a) 



for the lower bound, and 
Eu — ae 



f^l = % { [^HF - e(a - 1)] - ^[EuF - e{a - 1)]^ + Ae{E^ - Ehf)} , (7.3b) 

for the upper bound. Note that the solutions with the plus sign in front of the square root 
have been discarded, since they lead to an incorrect limit as e ^ 0. The energy derivatives 
of (17.31) read then 



d(3^ _ 1 
~dE ~2 



|[2Ei - e{a + 1) - Ehf] - \l[2Ei - t{a + 1) - Ehf]^ + 4e(E, - ^hf) 

^(T? \\^ 2Ei-ae-EuF \ . x 

+ {Ei - ae) < 1 . = } , (7.4a) 

/[2Ei - e{a + 1) - E^f? + MEi - ^hf) 
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for the lower bound, and 



^ |[£;hf - t{a - 1)] - y^l^HF - e{a - l)]^ + 4e(E„ - Ehf)} 



~{E-ae)—r (7.4b) 
^[EHF-e(a-iP + 4e(E„-i?HF) 



dE 2 



for the upper bound, which, when inserted into ( |6.10| ), provide an expression for the spec- 
troscopic factor, associated with the lower and the upper bound, respectively, in terms of 
(^o"),Ehf, aande. 

From eq. ( [7.3[ ) we note that inserting values of Ei (or E^) and of -Ehf one obtains a 
value of /3f and (3"^ for given values of a and e. These values of can then be inserted into 
eq. ( [5. 81 ) to obtain Eq and Eq. A solution of eq. ( |^ ) and ( |7.2| ) occurs when these upper 
and lower bounds for the mean field are equal. 

But how are Ei and E^ chosen? And how is the value of E'hf to be obtained? The first 
of these questions is answered as follows: 

Ei = S- W/2 (7.5a) 

Eu = S + W/2 , (7.5b) 

where 

^ = [-16 + 39.5(A;i.- 1.36)^]MeV (7.6) 

incorporates the present knowledge on the ground state energy of nuclear matter, quoted in 
ii), as extrapolated from finite nuclei. W is the fluctuation energy. 

The value of E'hf as a function of fci? is obtained from an assumed two-body interaction. 
This is described in Appendix 0. The resulting Hartree-Fock energy for a Fermi gas is 
shown in Fig. ^ and, for a reduced scale in kp, in Fig. |[ From Fig. ^ we see that for the 
assumed two-body interaction a minimum of the binding energy as a function of kp occurs 
at kp = 1.78 fm~^ with the corresponding energy per particle of —7.3 MeV. 

The potential can also be used to calculate the bare value of (3"^. This is discussed in 
Appendix |B[ Obtained values are given in Tables |, ^ and they range from 4.5 x 10^ to 
7.0 X 10^ MeVVnucleon. 

In discussing our findings we first observe that, for a given fluctuation energy W, a whole 
set of values for the parameters a and e exist such to satisfy the requirement i), namely 
Eq = Eq, at kp = 1.36 fm~^. They lie on the curves displayed in Fig. ^ Furthermore a and 
e should also be such to fulfill the constraint 

e2 = ae> 20 MeV , (7.7) 

which represents a fair estimate of the lower limit for the excitation energy of the 2p-2h 
states. This restricts the acceptable values for a and e to the domain to the right of the 
dotted line in Fig. ^ 

In this region a and e should be selected in such a way to comply with the requirements 
ii) and iii). This turned out to be possible and our results are shown in Fig. ^ and Tables |, 
fl|and|m[ 
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However there one sees that as one moves away from kp = 1.36 fm~^ the equahty 
between Eq and Eq is no longer exactly satisfied, although the two mean fields remain 
rather close in the range 0.9 < kp < 1.5 fm~^ providing the values of W are "moderate". 
Indeed, quantitatively, requirement i) can be reasonably satisfied in the above quoted range 
oi kp for W < 4: MeV only. For larger fluctuation energies Eq and Eq differ too much: an 
orientation on the size of the error around the mean field in thus obtained. 

Concerning the energy averaging parameter, e, is 1 MeV larger that W for all the cases 
listed in the tables. The values of the spectroscopic factors Si and Su are quite stable, while 
the values of the effective interaction jSf and (3^ grow rapidly with increasing W. Most 
importantly, their values differ by three orders of magnitude from the bare This is in 
part due to the renormalization (see eqs. (|3.13|) and (p.l4|) ) induced by the energy averaging. 



The random phase averaging also gives rise to a reduction in the magnitude of the residual 
interaction. Rough calculations indicate that these two effects are sufficient to produce the 
observed sharp reduction. 

VIII. CONCLUDING REMARKS 

The analysis presented in sections f|-fV| is based on two propositions. It is assumed that 
the mean field is the slowly varying component of the nuclear interaction, which can be 
obtained by taking an appropriate energy average. Secondly it is proposed that the matrix 
elements of the residual interaction are random, so that their average is zero. A formalism 
incorporating these ideas, borrowed from statistical reaction theory, was developed and 
explicit expressions for the mean field and the 'fiuctuation' away from the mean field was 
obtained. An expansion of the fiuctuation energy in terms of increasing excitation complexity 
leads, after averaging, to formulas for the corresponding contributions to fiuctuation energy. 

These formal discussion were followed in sect. and VI and [VII by a simplified version 



and by a schematic model. The most remarkable result is a sharp reduction the effective 
strength of the residual interaction. But in addition the reasonableness of our results is very 
encouraging. 

Of course much remains to be done. The quantitative connection with the underlying 
nuclear forces has not been exhibited. The evaluation of the matrix elements for finite 
nuclei was not carried out and a better understanding of the energy average needs to be 
achieved. We have so far only a schematic model. What is needed is a complete and thorough 
evaluation. But what has been indicated is that such an evaluation will be successful. 

Application to excited states is also indicated. In this case the smoothing function used 
in reaction theory can be used instead of the one of section |IV|. This would lead to a complex 
mean field and the excited state would have a width corresponding to the probability of the 
splitting of the state by the residual interaction. The width would then measure the extent 
of the splitting. 

APPENDIX A: THE BARE INTERACTION AND THE HF THEORY 



To implement the program outlined in Section VII we need a NN interaction to fix the 



HF energy of nuclear matter. For illustrative purposes we choose the following simple NN 
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interaction 



V(r) = QA 



9b 



(Al) 



r r 2 

which embodies a short-range repulsion in the first term and an intermediate range attraction 
in the second one (all the parameters are positive and we require //^ > /is). The latter is 
taken of Majorana type, hence the occurrence of the exchange operator 

1 + Ti ■ T2 1 



(Ti ■ 0-2 



2 2 

built out of the spin cr and isospin r operators. Thus ( |A1| ) contains the main features needed 
to account for the saturation of the nuclear forces, reflected in the existence of a minimum 
in the binding energy per particle (B.E./A) versus kp curve, but for the tensor force that 
here, for the sake of simplicity, is neglected. 

The HF energy for the interaction (^) is easily worked out and leads to the following 
expression for the binding energy per particle H] 



B.E. 

A 



3 h^kl Q 
5 2m ^ 2 



9a_ 




3kp 

4:71 



9A 



1 

+- 

TT 



gAl^A arctan 



iTiki 



3 + 



' 2kF ^ 



+ ■^gsfJ'B arctan 



log 1 + 



Aky 




(A3) 



+ 29Bf^B 



log 1 + 



where g = 2fc|,/37r^. 

The parameters characterizing the potential (|A1|) might be fixed, for example, by ac- 
counting for the "experimental" nuclear matter values previously quoted. One succeeds in 
doing so with the following choice 



(J'A 

9a 



3.43 fm"^ , 
2460 MeV fm . 



Hb = 1-63 fm ^ 
gB = 898 MeV fm . 



(A4a) 
(A4b) 



Worth noticing is that the range of the repulsion obtained with the fitting procedure is 
rather close to the one associated with the exchange of a a; meson (3.97 fm~^), whereas the 
range of the attraction turns out to be intermediate to the one arising from the exchange of 
a pion and of a a meson (0.71 fm~^ and 2.79 fm~^, respectively). 

Here, we rather prefer to choose the parameters in such a way to have too little binding 
energy at too large a density in the HF frame, in order to conform to a shortcoming common 
to many nuclear matter calculations, the purpose being to ascertain whether the present 
theory is capable to improve upon the HF results. Of course, there exists a variety of ways 
for reaching this scope: In view of the rather realistic values of the ranges ^a and fiB (see 
( |A4a|) ) we change the coupling constants. We thus take, as a rather extreme example. 



gA = 740 MeV fm , gB = 337 MeV fm , 



(A5) 



which, together with the values for /j^a and jj^s given in ( [A4a| ) , yields a minimum of —7.30 
MeV for the binding energy at kp = 1.78 fm~^, as it can be seen in Fig. |^, where the HF 
energy is displayed as a function of the density. 
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APPENDIX B: VACUUM ^ 2P 2H MATRIX ELEMENT 

In the language of second quantization the matrix element we have to calculate reads 
E E {k,,k2\V\k, + q,k2-q){k, + q,k2-q\V\k,,k2), (Bl) 

spin ki,k2<kp 
isospin |fci+<}|,|fc2— q|>A:F 



which, with standard manipulations, can be transformed into 

2 r dki dk2 
'~qJ (27r)3 (27r)3 



(3' = A- f ^J^dq Q{\k, + q\- M Q{\k,-q\- k^) 9 {k^ - k,) 

1 . . 1 



xe (fc. - k,) |32,,^-^^ + ^^sbj^^t^ - 24^-^-;^^ ;4T? 
/I 1 1 1 M 

-^'^QaQb — o— — ^ 2 , II 1 \ 1? + 



.lA + + 1^1 - ^2 + /i| + + 1^1 -k2 + q\ 



(the matrix elements in eq. ( |B1| ) are understood to be antisymmetrized) . The above 9- 
dimensional integral, by appropriate transformations, can be reduced to a combination of 2- 
and 4-dimensional integrals, which can be numerically evaluated, yielding the results quoted 
in Tables |-|I|. These values turn out to be helpful in performing a comparison with later 
findings. 

The direct contributions, namely the first three terms in the right hand side of (P^), are 
analytically evaluated. Considering, e. g., the third one, one gets 

A— / T^TTT^^^ Q{\k, + q\- M Q{\k2-q\- kp) e {kp - k,) e {kp - k^) 
f 1 1 1 

X 



'Akp {fii-I^B 



A— — (/^ A arctan yU A — /^B arctan yUs) (B3) 



4 /il - /i| 



/i'^(3 + yu|)^ arctan ( ) — /i^lS + /^a)^ arctan ( 



where /i^.b = fiA,B/'2kF. 

The other terms are obtained with obvious substitutions. 

The exchange contribution is somewhat more involved. As an example we consider (again 
the other pieces are obtained through obvious interchanges) 

2 f dki dk 

f^L = A— / TTT^TTr^dq {\k, + q| - kp) 9 (1^2 - q| - kp) 6 {kp - ki) 6 {kp - k^) 
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A 2 °° 

7^Zi-E(-l)'(2^ + l)/ d<l 
(27r^) n^g ^ Jo 



too 

/ dr re~'^'^'^ jo(qr) 
Jo 



(B4) 



X 



dkk^ji{kr) / dx Pi{x) e{q^ + k'^ + 2qkx - kp) 



where a partial wave expansion, in terms of the Legendre polynomials Pi and the spherical 
Bessel functions ji, has been performed. The above integral, for q > 2kF is easily calculated 
yielding 



Me: 



A 



Tl'^Q Jo 



dr — - — \ sm{kFr) — kpr cos{kFr 



2 r°° q 

dq ^ sin(gr) . (B5) 

2feF A«l + ? 



The remaining integrals are numerically evaluated. 



21 



TABLE I. In the second row the value of Eq as a function of kp is reported; in the third 
row one finds the bare value of as given by eq. (|B2|), while in the fourth and fifth rows the 
estimated renormalized values on the lower and upper bounds, respectively; in the last two rows, the 
corresponding values of the spectroscopic factor are displayed. The values in this table correspond 
toW = 2 MeV, e = 3 MeV and €2 = 22 MeV. 



kp (fm-i) 


1.2 


1.36 




1.5 


Eq (MeV) 


-15.3^-15.5 


-16.3 




-15.6H-15.2 




4.5-10^ 


5.2-10^ 




7.0-10^ 


3? fMeV^/nucleon) 


27.0 


26.3 




22.8 


Pi 


54.2 


49.0 




36.4 


Si 


0.32 


0.32 




0.30 




0.64 


0.62 




0.58 


TABLE IL As in table |: 


: The values in this table correspond to = 


3 MeV, 


e = 4 MeV and 


€2 = 20 MeV. 










kp (fm-i) 


1.2 


1.36 




1.5 


^0 (MeV) 


-15.5-1-15.9 


-16.6 




-15.9H-15.2 


P' 


4.5-10'^ 


5.2-10"^ 




7.0-10'' 


Pf (MeVVnucleon) 


36.1 


35.3 




30.7 


Pi 


79.6 


69.9 




49.8 


Si 


0.34 


0.33 




0.32 


Su 


0.69 


0.67 




0.63 


TABLE m. As in table p 


]: The values in this table correspond to W = 


4 MeV, 


e = 5 MeV and 


62 = 19 MeV. 










kp (fm-i) 


1.2 


1.36 




1.5 


Eo (MeV) 


-15.7^-16.4 


-16.8 




-16.1H-15.2 


P' 


4.5-10^ 


5.2-10'' 




7.0-10^ 


Pf (MeVVnucleon) 


45.7 


44.7 




39.2 


Pi 


108.6 


91.9 




62.7 


Si 


0.35 


0.34 




0.33 


Su 


0.73 


0.72 




0.66 
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FIG. 2. The binding energy per particle in the HF approximation (formula ( [A3| ) of the text) 
for the potential ( [Al| ) and the parameters given by ([A5|). 

12 I — I — '-I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — 




Also shown (dotted line) is the curve along which the average energy of the 2p-2h excitations is 
20 MeV. 
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-20 



_l I 1_ 



_l I 1_ 



1.2 



1.4 



1.6 




1.2 1.4 

FIG. 4. The binding energy per particle in the HF approximation (dot) and in the present 
approach: The sohd hnes represent the lower and upper bounds of fluctuations for the energy E, 
whereas the dashed lines give the corresponding mean field energies, (a): W = 2 MeV, e = 3 MeV, 
a = 7.3, £2 = 22 MeV; (b): W = 2, MeV, e = 4 MeV, a = 4.9, 63 = 20 MeV; (c): W = A MeV, 
e = 5 MeV, a = 3.8, €2 = 19 MeV. 



24 



REFERENCES 



[1] H. Feshbach, "Theoretical Nuclear Physics: Nuclear Reactions", J. Wiley & Sons, New 
York (1992). 

[2] O. Benhar, A. Fabrocini and S. Fantoni Phys. Rev. C41 (1989), R24. 
[3] H. Feshbach, Phys. Rep. 264 (1996), 137. 

[4] M. Kaway, A.K. Kerman and K.W. McVoy, Ann. Phys. (N.Y.) 75 (1975), 156. 

[5] H. Feshbach, A. K. Kerman, and S. Koonin, Ann. Phys. (N.Y.) 125 (1980), 429. 

[6] E. K. U. Gross, E. Rung, and O. Heinonen, "Many-particles theory", Adam Hilger, 

Bristol (1991) p.59. 
[7] H. A. Bethe, Ann. Rev. Nucl. Sci. 21 (1971), 93. 

[8] W. M. Alberico, M. Bregola, R. Cenni, and A. Mohnari, Nuovo Cimento 58A (1980), 
31. 



25 



